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A gamma Stirling engine is studied in this paper. A special care was accorded to the instrumentation of 
this engine and especially the instrumentation of the regenerator. A preliminarily set of experimental 
measurement reveals a difference of temperature between both regenerator sides. A second set of exper¬ 
iments was proposed to detect the influence of this phenomenon on Stirling engine performances. The 
asymmetry of heat transfer inside the Stirling engine regenerator’s is one of the important phenomenons 
which consume a part of the produced energy. Two experiments are made to find out the causes of this 
asymmetry. In order to know the influence of the different operation parameters on this new phenome¬ 
non the experimental design method is adopted. The experimental design is an alternative to identify the 
parameters sets allowing optimal Stirling engine performances. A central composite rotatable design was 
adopted for minimizing the asymmetry of temperature between both regenerator sides and maximizes 
the engine brake power. The selected four independent parameters are: heating temperature (300 °C- 
500 °C), initial filling pressure (3 bar-8 bar), cooling water flow rate (0.2 l/m-3 1/min) and operation time 
(4-20 min after study regime). The four adopted factors are experimentally varied. The results show that 
the heating temperature is the most significant factor for the studied phenomenon. The major damages 
caused by this phenomenon will be presented too. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

The optimization of Stirling engine performances is one of the 
priorities of the related researchers in the past several years. So, 
many experimental and analytic studies were done to understand 
phenomena of heat transfer [1], irreversibility [2], imperfect regen¬ 
eration and flow friction [3] in Stirling engine. The performance of 
such engines is very sensible to its geometrical [4-7 and operation 
conditions [6-10 . 

The temperature distribution in Stirling engine is different from 
a compartment to another. The most important variations are 
observed in the regenerator (porous media, complicated heat 
exchange, heated and cooled at the same time). Temperature 
distribution is a determinant parameter to the regenerator effec¬ 
tiveness. The regenerator is also the seat of an important part of 
thermal losses recorded in a Stirling engine [11 . The most consid¬ 
ered losses are lost by pressure drop, power lost by external and 
internal conduction and the regenerator ineffectiveness. Tlili 
et al. [11] showed that the power loss observed in the Stirling 
regenerator are about 86% of the total engine losses. 
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Popescu et al. [2] propose a finite time thermodynamic optimi¬ 
zation of a Stirling engine based on endo and exo-irreversible cycle. 
They conclude that the internal irreversibilities are due to the 
internal conductance of the plant and to the non-adiabatic regen¬ 
erator (drops of 30% of engine performances). The external irrever¬ 
sibilities are due to the linear interactions between the heat 
sources and the working fluids at finite temperature gaps. 

The regenerator parameters have been the subject of several 
experimental and analytical investigations. The regenerator 
effectiveness was one of the decision variables used in the Multi¬ 
objective evolutionary algorithms developed by Ahmadi et al. 
[12]. They proposed an optimal Pareto frontier in objectives spaces. 
Chen et al. [13] built a prototype helium-charged twin-power pis¬ 
ton c-type Stirling engine. Some of its geometrical and operational 
parameters were investigated by a numerical model. They found 
that regeneration effectiveness had the most prominent effect on 
efficiency, while engine speed had the greatest effect on the engine 
power within the range of the engine speed. Cheng et al. [14] devel¬ 
oped and tested a beta-type 300-W Stirling engine. The engine 
parameters are insured in a non-ideal adiabatic model to predict 
its performances. Using a 120 wire mesh regenerator, the shaft 
power of the engine reaches 390 W at 1400 rpm with 1.21-kW 
input heat transfer rate (32.2% thermal efficiency). Kato and baba 
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Nomenclature 



bo 

constant coefficient 

W 

diagonal matrix 

bj 

linear coefficients 

X 

level matrix of independent variable 

b„ 

quadratic coefficients 

Xi 

coded values of the factor or regression coefficient 

by 

second order interaction coefficients 

Xcp 

independent variable real value at the center point 

m 

cooling water flowrate 

X* 

average value of X* 

N 

total number of experiments 

A Xi 

the step change of the real value corresponding to a var- 

Pi 

initial filling pressure 


iation of a unit for the dimensionless value of the vari- 

a 

degrees of freedom 


able i 

T 

temperature 

Y 

vector of observation 

T h 

heating temperature 

Ymod 

predicted response 

Top 

operation time 

y* 

average value of Y 

a r exp 

experimental difference of temperature between both 




regenerator sides 

Greeks symbols 

A T pr 

difference of temperature between both regenerator 

a 

critical value of Pareto chart 


side obtained from the predicted model 

p 

regression vector coefficients 

u, 

experimental response 

fi 

mean value of the x variables 

urn 

value of Ui in the middle 

o 

standard deviation 

A Ui 

maximum variation of U { 




[15] proposed an experimental method for regenerator evaluation. 
They concluded that the efficiency of the regenerator depends on 
the difference of temperature of the working fluid on the hot side 
and on the cold side and the pressure fluctuation. Chen et al. [16] 
noticed that the operating temperature ratio dependent on the 
thermal resistance of heating and cooling spaces, the regenerator 
effectiveness, the working fluid mass and the engine speed. Glu¬ 
shenkov et al. [17] proposed a thermodynamic analysis based on 
a heated regenerator. The developed thermodynamic model incor¬ 
porates non-ideal features of the cycle, such as specific regenerator 
efficiency, dead volumes and other geometrical parameters. The 
model shows that the energy efficiency is highly sensitive to regen¬ 
erator performance. Ahmadi et al. [18], optimize a solar-powered 
high temperature differential Stirling engine with multiple criteria. 
The volumetric ratio and the regenerator effectiveness effects on 
the final selected solution were studied and it was found that a final 
optimal solution with higher values of objectives (thermal effi¬ 
ciency and output power) could be achieved if these parameters 
were increased. Kato et al. [19] presented an empirical estimation 
of Stirling engine regenerator efficiency. The regenerator efficiency 
is calculated as a function of the temperature fluctuation in the 
regenerator and the temperature fluctuation of cold side working 
fluid. They showed that these two parameters are positively corre¬ 
lated with the regenerator loss. Cheng and Yu [20] proposed a more 
realist model to study beta-type Stirling engines. Their model took 
into consideration the temperature difference between the working 
gas and the external heat source, pressure drop in the regenerator, 
and regenerator effectiveness. Formosa and Despesse [21 ] included 
dead volume, imperfect regenerator, and external and internal ther¬ 
mal transfers in their analytical model that could be used to analyze 
Stirling. Chen et al. [22] experimentally studied the effects of sev¬ 
eral regenerator parameters on the overall performance of Stirling 
engine. They incorporated a moving regenerator in their engine. 
The parameters investigated include regenerator matrix material, 
matrices arrangement, matrix wire diameter, and fill factor. Ander¬ 
son et al. [23] concluded that the convective heat transfer is the 
dominate heat transfer mechanism between regenerator matrices 
and working gas. Regenerator heat transfer rate is a very crucial fac¬ 
tor to regenerator effectiveness especially at high engine speed 
when the time for heat transfer is very short in each cycle. There¬ 
fore, a regenerator with low heat transfer rate will perform poorly 
at high engine speed even though it has large heat capacity. Gheith 


et al. [24 experimented a gamma type Stirling engine with differ¬ 
ent regenerator materials (Cooper, Aluminum, Stainless Steel and 
Monel 400). The obtained experimental results provide guidance 
to Stirling engine enhancement and selection of the appropriate 
regenerator constituting material. The Monel 400 and the stainless 
steel regenerators present an acceptable thermal efficiency and do 
not oxidize. Tavakolpour et al. [25] studied numerically a LTD (Low 
Temperature Difference) Stirling engine. They found that the 
energy losses in the regenerator are due to temperature oscillation 
and heat conduction caused by the alternating gas flow and steep 
temperature gradient throughout the regenerator and found that 
the engine’s thermal efficiency resulted from the regenerator 
efficiency of 1.0 at ideal conditions the is 0.069 whereas for the 
regenerator efficiency of zero, is about 0.0122. The use of an effi¬ 
cient regenerator increases six times the engine thermal efficiency. 

The experimentation of Stirling engine is always made by clas¬ 
sical way (univariate study) [25-30 . The investigation by such 
method requires a higher number of experiments for exploring 
the whole experimental domain. It is time consuming and does 
not depict the complete effect of each parameter in the response 
(Stirling engine performances). 

To reduce these disadvantages during experimental studies, 
especially for study involving several variables factors, the experi¬ 
mental design methodology seems to be adequate to avoid such 
inconvenient. 

The experimental design methodology is not a new technique it 
was developed by Fisher in 1952 [31] and applied first in 
agronomy. -This method presents several advantages over the 
conventional one: 

- Decrease the number of test performed to scan the entire field 
of study. 

- Possibility to study a large number of factors at a time. 

- Detect single or double interactions between the studied 
factors, [32]. 

- Modeling the studied responses. 

- Determining optimal results [31 . 

The aim of this paper is to study a new loss which dissipates a 
part of Stirling engine performances. Detailed descriptions of the 
experimental set-up as well as its corresponding metrology are 
presented. 
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The first section presents a global characterization of the 
temperature distribution on the regenerator. A zoom on both 
regenerator sides shows an important gap of temperature. By per¬ 
forming additional experimental measurements it was possible to 
determine the causes and the consequences of this phenomenon 
on the engine performances. 

A new way of optimization of Stirling engine performances is 
proposed in the second section of this paper with consideration 
of the new observed phenomenon. The experimental design meth¬ 
odology was adopted to detect the optimal parameters minimizing 
the anomaly of temperature distribution in the regenerator and 
consequently maximizing the performances of the Stirling engine. 
For operation factors (heating temperature, initial filling pressure 
and cooling water flow rates) are varied and their corresponding 
influence on the studied responses are detected and discussed. 

2. Experimental set up 

2.2. Gamma Stirling engine 

The experimental setup is presented in Fig. 1. It is a gamma 
Stirling engine, working with a maximal compressed air of 
lObar. The maximal rotation speed is around 600 rpm and 
provides 500 W of mechanical shaft power. In addition to the 
expansion and compression spaces, the engine includes three heat 
exchangers: the regenerator, the cooler, and the heater. The regen¬ 
erator (Fig. 1) is a porous medium (i.e., made of copper with 90% 
porosity). The cooler acts as a cold source formed by an open water 
circuit, and the heater consists of 20 tubes in order to increase the 
exchange surface between the working gas and the hot source 


Table 1 

Technical specifications of the Stirling engine. 


Geometrical parameters 

Compression space, “c” 


Regenerator, “r” 


Diameter 

65 mm 

Diameter 

52 mm 

Height 

170 mm 

Height 

153 mm 

Expansion space, “e” 


Porosity 

0.9 

Diameter 

65 mm 

Heater, “ft” 


Height 

170 mm 

Diameter 

12 mm 

Cooler, “/<” 


Height 

650 mm 

Diameter 

12 mm 

Crank radius 

66 mm 

Height 

900 mm 

Compression ratio 

3.5 

Functioning parameters 
Nominal rotation 

600 rpm 

Working fluid 

Air 

speed 




Hot source 

500 °C 

Cold source 

15 °C 

temperature 


temperature 


Porosity 

0.9 

Regenerator material 

Stainless 

steel 


(electrical resistance). The operating conditions and the geometri¬ 
cal characteristics of the studied Stirling engine are recapitulated 

in Table 1. 

2.2. Instrumentation 

To characterize the Stirling engine, we performed the experi¬ 
ment using different parameters such as the pressure, the temper¬ 
atures at different positions in the engine, and the cooling water 
flowrate. Several type-K thermocouples with a diameter of 0.5, 



Fig. 1. The Stirling engine set-up: 1 - expansion space, 2 - drive system, 3 - cooler, 4 - regenerator, 5 - compression space, 6 - Heater, 7 - heating system, 8 - dissipation 
system. 
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4 Skinned 
thermocouples 



4 Skinned 
thermocouples 


Fig. 2. The regenerator instrumentation. 


0.25, and 0.0254 mm were used to measure the temperature at dif¬ 
ferent locations of the engine. Because the interest was to measure 
the instantaneous temperatures, proper care has been taken 
regarding the size and the thermal inertia (which should be very 
low) of the thermocouples. Two thermocouples were mounted 
upstream of the expansion and compression spaces, respectively. 
Other two thermocouples were used to measure the inlet and 
outlet temperatures of the water supplying to the cooler. Eight 
thermocouples are implemented to measure the temperature of 
the working gas passing through the symmetrical side of the 
regenerator matrix. The thermocouples were skinned up to 1 mm 
in the regenerator matrix without touching the material (Fig. 2). 
The brake power produced by the Stirling engine was transmitted 
by a transmission belt to an alternator. The alternator was 
mounted on a balance plate that has been used to determine the 
torque transmitted from the transmission belt to the alternator. 
A flow-meter localized upstream of the entrance of the engine in 
the water-cooling circuit was used to determine the flow rate of 
the cooler. A crank angle transducer was connected to the engine 
shaft. Pressure transducers were installed in the compression 
space and after the expansion space. Their frequency response 
was chosen for a wide range to ensure good instantaneous mea¬ 
surements of the pressures inside the working spaces. 

Acquisition data were recorded using a fast acquisitions system 
(ADwin Pro-II, Jager Computergesteuerte Messtechnik GmbH, Lors- 
ch, Germany). The ADwin Pro-II is an external processing system 
with modular expansion options. Our system includes eight mod¬ 
ules in parallels. The first two modules were used to the acquisition 
of the thermocouple signals. Each module is formed of eight mul¬ 
tiplexes input with 18 bits of precision. Hence, the maximum 
uncertainties on the temperature measurements do not exceed 
0.002 °C. The transducers of the pressure in each of the two space 
rooms, the force and the flowrate, are linked to two analog input 
modules. Each one of these modules has eight fast analog-to- 
digital converter of 18 bits and eight differential inputs. The digital 
input/output module provides 32 programmable digital input and 
output channels at transistor-transistor logic levels. All these mod¬ 
ules are controlled by the processor module. It has an internal 
memory of 768 KB and has an external trigger input that enables 
the processor module to recognize an external signal. A crank angle 
transducer is located in the end of the Stirling engine shaft. On 
every crankshaft degree, this crank angle transducer sends a tran¬ 
sistor-transistor logic external signal to the processor module, 
which allows the acquisition of all signals by all used modules. 
For every cycle, 360 acquisitions are done, which correspond to 
360° of the crankshaft. For every run, 100 cycles are recorded. 
The acquisition duration depends on the Stirling engine revolu¬ 
tions per minute (rpm) for the entire 36,000 acquired measure¬ 
ment points and for each signal. 

3. Global study of the regenerator distribution of temperature 

A Stirling engine regenerator is crossed twice by the working 
fluid during each Stirling cycle. Assume that the direction of the 


working fluid flow is from the cooler to the heater through the 
regenerator (porous media). During the first phase of the cycle, 
the regenerator porous matrixes discharge its stored heat to the 
working fluid (preheating the working fluid). In the opposite direc¬ 
tion of flow (the working fluid flowing from the heater to the 
cooler), the porous matrix is heated by the working fluid (pre-cool¬ 
ing of the working fluid). 

Thanks to the adequate thermocouples (very small diameters) 
placed on both regenerator sides, the instantaneous temperatures 
of the working fluid at different positions of the regenerator are 
recorded and represented in Fig. 3. It can be seen that all detected 
temperatures present periodic evolution with different mean 
value and amplitude. The temperature distribution in both sides 
of the regenerator matrix is not identical. Side 1 of the regenerator 
matrix is warmer compared to side 2. This observation is con¬ 
firmed by the axial velocity recorded on both regenerator sides. 
This parameter was calculated using the shift phase between T R \ 
and T R4 for side 1 and between T R5 and T R8 for side 2. The axial 
velocity of the working gas in the regenerator side 1 is about 
0.23 m/s. However, it is about 0.46 m/s in the side 2. The axial 
velocity of the working fluid in the side 2 of the regenerator is 
twice than that of side 1. Hence, the heat exchange duration is 
more important in side 2. The slowest working fluid (side 2), 
spends more time in the regenerator and therefore has more time 
to interact with the porous matrix and absorb more calories. 
It is clear that the heat transfer in both regenerator sides is 
asymmetrical. 

The experimental working gas velocities in both regenerator 
sides were calculated using the shift phase between T R i and T R4 
for side 1 and between T R5 and T R8 for side 2. The axial velocity 
of the working gas in the regenerator side 1 is about 0.23 m/s. 
However, it is about 0.46 m/s in the side 2. The axial velocity of 
the working fluid in the side 2 of the regenerator is twice than that 
of side 1. Hence, the heat exchange duration is more important in 
side 2. In fact, the working fluid spends more time in side 1 absorb¬ 
ing more energy. It is clear that the heat transfer in both regener¬ 
ator sides is asymmetrical. 

The heating temperature and the initial filling pressure were 
varied and the difference of temperature between two thermocou¬ 
ples placed on opposed regenerator side but at the same level are 
recorded. The difference of temperature between both regenerator 
sides increase with the heating temperature. It has an average 
value of 37 °C for a heating temperature of 500 °C (Fig. 4). The ini¬ 
tial filling pressure does not have an important influence on the 
studied phenomenon. For an increase of 5 bar, the difference of 
temperature increases of only 2 °C for all applied heating 
temperature. 

It is a new lost, which dissipate an important part of thermal 
energy exchanged in the regenerator. This lost is not considered 
in the theoretical models encountered in the literature. It depends 
mainly of the compartment dispositions and design. In the follow¬ 
ing section the causes and the influence of this lost on the engine 
performances will be discussed. In addition of internal conduction 
loss and lost by imperfect regeneration, the lost by dissymmetry 
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Fig. 3. The working fluid temperature evolution at different position of the regenerator. (Thermocouples T R i to T R4 in side 2 & Thermocouples T R5 to T R8 in side 1). 
Experimental conditions: T h = 500 °C, P, = 8 bar, rh = 8.16 1/min. 


40 
38 
36 
34 
32 

o 

“ 30 
< 

28 
26 
24 
22 
20 

Fig. 4. Evolution of the mean difference of temperature between both regenerator 
sides function of the initial filling pressure and for different heating temperatures. 



regenerator of temperature in the can be added t theoretical model 
in same case of Stirling engine. 


cooling water is at its minimum temperature, coming to the oppo¬ 
site side of the regenerator, the water becomes hotter. 

The previous conclusion is confirmed by a first experience 
(Fig. 6A). All regenerator temperatures are records when the 
engine is stopped, and heated at 100 °C. Initially, all thermocouples 
record nearly the same temperature, a difference is gradually 
recorded between all thermocouples and especially between both 
regenerators’ sides. After 5300s, the cooling water flow rate is 
opened; immediately all recorded temperatures decreases. The 
temperatures T R i and T R5 (the nearest to the cooler) are different, 
the working fluid side 1 (T^i thermocouple) of the regenerator is 
lower than that of the fluid at the same level but in the opposite 
side 2 (thermocouple T R5 ). 

The Gamma Stirling engine was heated at a temperature of 
300 °C, reaching this temperature the engine was started. All 
recorded temperatures are represented in Fig. 6B. Initially, the 
thermocouples present almost the same temperature of the work¬ 
ing fluid on both regenerator sides. Then (engine starts), the tem¬ 
peratures have increased over time and eventually stabilize when 
the engine has reached its quasi-steady operation. Temperatures 
delivered by thermocouples T R i and T R5 , increased with the engine 
started, but the gap between them has also increased. Whit the 
engine start this difference increases from 1 °C to 20 °C on the 
quasi-study regime. 


4. Causes of the asymmetry of temperature between both 
regenerator sides 

It is clearly showing in Fig. 5 that the engine design is asymmet¬ 
rical. The circulation of the working fluid from the compression 
working space to the cooler is made on the left side of the engine 
(side 1). The entrance of the cooling water (dark blue) is made on 
the right side of the Stirling engine (side 2). The side 1 (right side) 
of the regenerator is more cooled by the water flow rate (heat 
transfer by conduction through the wall separating the cooler from 
the working fluid). Indeed at its input to the Stirling engine, the 


5. Experimental design methodology 

The experimental design methodology was used to study the 
effect of four operation parameters on the asymmetry of tempera¬ 
ture between both regenerators’ sides. The heating temperature, 
the charge pressure, the cooling water flow rate and the operation 
time were chosen as independent factors. 

The Yate notation (Fig. 7) allows a simple representation of 
combinations of levels and allows easy calculation of factors effects 
and interactions. Levels are represented by + and - signs symboliz¬ 
ing the high and the low value of a factor. The relationship between 
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Working fluid input 


Cooler 


Cooling water 
output 



Series of 

thermocouples on -<( 
regenerator side 2 


Cooling water input 


Series of 
thermocouples on 
regenerator side 1 


Expansion 

space 


Fig. 5. Part of the cut of the studied Stirling engine. 


A 



B 



Fig. 6. (A and B) Experiments to determine the existence of the difference of temperature between both regenerator sides. 


the coded factors and the encoded ones are calculated by the fol¬ 
lowing formulas: 


,/ _ X\ X C p 

1 ~ A x t 




X-imax 

2 



In order to detect the influence of each studied factor, two 
experimental designs were conducted: the first was a Factorial 
Centered Design (FCD) for determining the most influencing fac¬ 
tors and the second was a Central Composite Rotatable Design 
(CCRD) for determining the response surfaces model for the more 
influencing factors and then obtaining the optimal functioning 
conditions. 
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Fig. 7. Yate notation. 


The multiple regression analysis was used to obtain a qua¬ 
dratic response surface. The empirical model linking the difference 
of temperature between both regenerator sides to all studied 
factors Xi , their interactions x, Xj and the quadratic terms x 2 is 
writing as: 

4 

^mod — ^0 H - ^ ^ jjXj "T ^ ^ bjkXjX-k (3) 

j= 1 0<j</<<4 

Eq. (3) can be expressed in the matricial form as follows: 

T = XB + s (4) 


The sum of squared error (SSE) is defined by: 

SSE = (Y - XB) t (Y - XB) (5) 

Taking derivative of SSE with respect to B gives the minimum of 
SSE from the normal equation: 

X t XB = X t Y (6) 

So, the estimate, b , of the elements of B is given by: 

b= (x T -Xy\ T Y (7) 

The total sum of square can be partitioned into two parts, the 
sum of square due to regression (SSR) and the sum of square unac¬ 
counted for by the fitted model (SSE): 

SSR = (Y — Xb) [Y - Xb) T (8) 

With 

(9) 

With consideration of Eq. (4), the total sum of squares is written 
as: 

SST = SSR + SSE 

= l (Y ~ Xb) (Y-Xb) T + ± (Xb - Y)(Xb - Y) t (10) 


Table 2 

Coded and actual levels of factors considered for design. 


Factors 

Levels 





Lowest (-2) 

Lower (-1) 

Center(0) 

High (+1) 

Highest (+2) 

(A) Cooling flow rates (1/min) 

0.26 

2.25 

4.21 

6.17 

8.16 

(B) Charge pressure (bar) 

3 

4.25 

5.5 

6.75 

8 

(C) Heating temperature (°C) 

300 

350 

400 

450 

500 

(D) Operating time (s) 

4 

8 

12 

16 

20 


Table 3 

Design and the yield response. 

Run 

Cooling flow rates (1/min) Charge pressure (bar) Heating temperature (°C) Operating time (s) Ar exp (°C) 

Factorial points 






1 

2.25 

4.25 

350 

8 

32.28 

2 

6.17 

4.25 

350 

8 

28.10 

3 

2.25 

6.75 

350 

8 

29.29 

4 

6.17 

6.75 

350 

8 

28.34 

5 

2.25 

4.25 

450 

8 

42.90 

6 

6.17 

4.25 

450 

8 

37.70 

7 

2.25 

6.75 

450 

8 

41.90 

8 

6.17 

6.75 

450 

8 

36.68 

9 

2.25 

4.25 

350 

16 

34.28 

10 

6.17 

4.25 

350 

16 

28.10 

11 

2.25 

6.75 

350 

16 

30.29 

12 

6.17 

6.75 

350 

16 

28.34 

13 

2.25 

4.25 

450 

16 

43.90 

14 

6.17 

4.25 

450 

16 

37.70 

15 

2.25 

6.75 

450 

16 

40.03 

16 

6.17 

6.75 

450 

16 

36.68 

Centre points 






17 

4.21 

5.5 

400 

12 

33.97 

18 

4.21 

5.5 

400 

12 

33.01 

19 

4.21 

5.5 

400 

12 

34.53 

Axial points 






20 

0.26 

5.5 

400 

12 

28.50 

21 

8.16 

5.5 

400 

12 

26.14 

22 

4.21 

3 

400 

12 

30.40 

23 

4.21 

8 

400 

12 

29.87 

24 

4.21 

5.5 

300 

12 

24.70 

25 

4.21 

5.5 

500 

12 

47.84 

26 

4.21 

5.5 

400 

4 

32.58 

27 

4.21 

5.5 

400 

20 

35.88 
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And then the coefficient of multiple determinations R 2 and the 
adjusted coefficient of multiple determinations R ad 2 can be 
determined as follows: 

CCD 

R 2 = — = 1 - SSE/SST (11) 



SSE/(n — k — 1) 
SST/(n - 1) 



The experimental design must choose a set of n combinations of 
the studied factors. The set of pilot experiment must be different 
from each other and must belong to a pre-defined domain to 
explore. The combinations are chosen to efficiently investigate 
the relationship between the design factors and the response. 
The experimental plan generated using the Minitab 16 Software, 
along with the results is depicted in Table 3. 

5.2. Determination of the significant parameters on the studied 
response 


In order to determine the exact term (influent factors, interac¬ 
tion and quadratic terms) of the predicting model, the following 
approach is used: 

The Pareto chart allows detecting the factors and the interac¬ 
tions which are most important to the studied response [33]. In 
this chart, the length of each bar is proportional to the absolute 
value of the estimated effect [34]. Effects that cross the vertical line 
(y = a ) are statistically significant at 95% confidence level. The 
value of a is calculated from the t- test [35]: 


a = 


(X* - M) 
o ■ yh 



The single effect of all studied parameters as well as their inter¬ 
actions (cross effects) up through order 3 can be discussed from the 
Pareto chart illustrated by Fig. 8. The numerical estimates of the 
effects indicate that the heating temperature (C) has the largest 
effect on the difference of temperature between both regenerator 
sides. The charge pressure (B) and the cooling water flow rate 
has a significant effect too. The operation time an all considered 
interactions do not have influence on the studied response. 
Certainly, the Pareto chart can detect the meaningful factors and 
interactions, but it cannot determine the significance of quadratic 
interactions (A 2 , B 2 , C 2 and D 2 ). 


2,45 



Standardized effects 

Fig. 8. Pareto chart. A: cooling water flowrate, B: initial pressure, C: heating 
temperature and D: operation time. 


Table 4 

Critical probability test. 


Variables 

Critical probability P 

Decision significance with an error of 5% ? 

bo 

0.000 

(Hi) Yes 

A 

0.000 

(HA Yes 

B 

0.000 

(HA Yes 

C 

0.000 

(HA Yes 

D 

0.137 

(H 3 ) Yes 

A 2 

0.011 

(H 2 ) Unknown 

B 2 

0.000 

(HA Yes 

C 2 

0.062 

(H 3 ) No 

D 2 

0.005 

(H 3 ) Unknown 

AB 

0.786 

(H 3 ) No 

AC 

0.689 

(H 3 ) No 

AD 

0.997 

(H 3 ) No 

BC 

0.737 

(H 3 ) No 

BD 

0.796 

(H 3 ) No 

CD 

0.699 

(H 3 ) No 


Table 5 

Student test for factors and interactions. 


Variables 

T-Student 

T-Student critique 

Decision 

A 2 

-2.894 

2.131 

Significant 

D 2 

-3.325 


Significant 


To determine the influence of quadratic interactions (A 2 , B 2 , C 2 
and D 2 ) on the studied phenomenon, the critical probability test 
is used. It is called degree of significance too. It indicates the per¬ 
centage of risk take when choosing a sample. Each sample is com¬ 
pared to the error value (5%). To make a decision on the 
significance of a factor, three hypotheses are proposed. 

Hi. P < 0.01 the factor is significant for the response. 

H 2 . 0.01 < P < 0.05 Unknown result, must use the Student and 
the Fisher tests. 


H 3 . P > 0.05 the factor is not significant for the response. 

Table 4 depicted all variables and their corresponding decisions. 
We found that the test of the critical probability is not 
sufficient to judge the significance of the variables A 2 and D 2 . 
So, we’ll go to the Student test to assess the significance of these 
two variables. 



> t(a, n 


Q 




The following table ( Table 5) summarizes the values of Student 
variables, and the decision for each of them. The t critical is deter¬ 
mined from statistical table found in literature. 

Data were analyzed using multiple regression analysis. So, with 
consideration of the significant variables for the response: asym¬ 
metry of temperature between both regenerator sides, the empir¬ 
ical model as function of coded variables can be written as: 


Y mod = 32.1476 - 1.325871 + 1.0817B + 4.595C - 0.6286A 2 

- 1.0161B 2 -0.7224D 2 (15) 

With consideration of Fable 2 and the Eq. (2), the prediction model 
function of uncoded moded is writing as follows: 

Y* mod = -38.4467 + 0.8276m + 7.998P/ + 0.0919 T H 

- 0.2227m 2 - 0.6503 Pf - 0.045ir 2 p (16) 
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5.2. Regression test 

The coefficient of multiple determinations R 2 and the adjusted 
coefficient of multiple determinations R ad 2 are respectively of 
0.9721 and 0.9442. Both coefficients are statistically significant 
for the studied response. It suggests that the estimated regression 
equations for the Case Studied fit the data well. Predicted and 
experimental results are represented in Fig. 9. It can be seen that 
all points cluster around the line y = x. This change gives an idea 
of the obtained predicted model and shows that it is in good agree¬ 
ment with the experimental values. 



Fig. 9. Normal probability plot of residuals for dissymmetry of temperature 
between both regenerator sides. 


5.3. Iso-surface response 

A technique used to help visualize the shape of the three- 
dimensional response surface is to plot the contours of the iso¬ 
surface response. In a plot contour, lines or curves of constant 
response values are drawn on a graph or plane whose coordinate 
axes represent the levels of the factors. The contour plot can be 
used to visualize the response surface behavior function of two 
factors at a time. Four factors are studied, so for every plot two 
factors will be hold at a constant level and plotting the other two 
factors. The iso-response surfaces for the factors (the heating tem¬ 
perature, cooling flow and initial charging pressure) influencing 
the asymmetry of temperature in the regenerator are plotted. They 
have been plotted in pairs. Fig. 10A. indicates that how variables, 
heating temperature and initial filling pressure, are related to the 
studied response while the other factors, cooling water flow rates 
and operation time are held constant at the highest level (+2). 
The difference of temperature between both regenerator sides is 
at its highest value (greater than 30 °C) at the darkest region of 
the graph. For low heating temperatures this difference of temper¬ 
ature is low (<20 °C) for all variations of the initial charge pressure. 

For a heating temperature of 350 °C (level-1), we keep the 
same level of gap of temperature between both regenerator side, 
but beyond this level this difference increases significantly, reach¬ 
ing a maximum about 30 °C for a heating temperature of 500 °C 
(level +2). For low cooling water flow rate (Fig. 10B), the difference 
of temperature between both regenerator sides increases of 
approximately 35 °C for a heating temperature of 500 °C (level + 2). 
The increase of cooling water flowrate decreases the difference of 
temperature between both regenerator sides. It can be noted in 



-l 


-2 -\ - 1 - \ - 1 - 

- 2-10 12 

Flow rate 


Fig. 10. Contour plot of the difference of temperature between both regenerator sides. 
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Table 6 

Operation parameters for minimal difference of temperature. 


Factors 

Uncoded value 

Heating temperature 

300 °C 

Initial filling pressure 

3 bar 

Cooling water flow rate 

8.16 1/min 

Operation time 

20 min 


Fig. 9C that regardless of the cooling water flowrate, the difference 
of temperature between both regenerator sides is very high. The 
increased of cooling water flow rates from the lowest level (-2) 
to its highest level (+2), causes the decrease of the temperature dif¬ 
ference of 12 °C. 

From the predicted model (Eq. (16)), it is possible to determine 
the optimum parameters minimizing the asymmetry of heat trans¬ 
fer between both regenerator sides. It was found that a minimal 
difference of temperature can be obtained for the operation param¬ 
eters depicted in Table 6. 

The asymmetry of temperature in a given section of the regen¬ 
erator is the consequence of air recirculation zones at the outlet or 
at the inlet of the porous matrix. It causes overlap (friction) 
between the molecules of air at different temperatures, and there¬ 
fore at different speeds. Consequently, a frictional drag force 
appears causing the dissipation of mechanical energy on heat. A 
difference of density between the particles of warm and cold work¬ 
ing fluid causes a pressure drop in the Stirling engine. The pressure 
drop causes deterioration of heat exchange in the engine, and 
hence a degradation of performance of the engine. Therefore, 
stationary vortexes appear at the entrance or the outlet of the 
regenerator. These vortexes cause the slowdown of the flow in 
some areas inside the porous matrix. 

6. Multi-objectif optimization 

The influence of the studied parameters on the engine perfor¬ 
mances with consideration of the dissymmetry of temperature 
between both regenerator sides is the subject of a multi-objectif 
optimization. The goal of this optimization is to minimize the dis¬ 
symmetry of temperature and to maximize the engine brake 
power. The Optimization Criteria are depicted in Fable 7. 

The superposed surfaces are a superposition of the iso-surfaces 
of both studied responses. The common part (white zone) visualize 



Initial filling pressure 


Fig. 11. Superposed surface for both studied responses function of initial filling 
pressure and heating temperature. 



Cooling water flow rate 


Fig. 12. Superposed surface for both studied responses function of heating 
temperature and cooling water flow rates. 



Cooling water flow rate 


Fig. 13. Superposed surface for both studied responses function of initial filling 
pressure and cooling water flow rates. 


the maximal obtained power and the minimal obtained difference 
of temperature with the variations of a couple of operation param¬ 
eters (heating temperature and initial filling pressure (Fig. 11), 
heating temperature and cooling water flow rates (Fig 12) and ini¬ 
tial filling pressure and cooling water flow rates (Fig 13). 

It can be seen (Fig. 11) that a maximal initial filling pressure 
associated to an average value of the heating temperature gives a 
brake power between 200 W and 290 W and a difference of tem¬ 
perature between both regenerator sides between 25 °C and 35 °C. 

It can be seen in Fig. 12, that the cooling water flow rates and 
the heating temperature has a significant influence on both studied 
responses. A maximal flow rates associated to mean temperature 
(400 °C) gives the prompted results (choosing in iable 6). 

A maximal initial filling pressure associated to a maximal cool¬ 
ing water flow rates gives a maximal brake power. This two oper¬ 
ation parameters are not very influent for the difference of 
temperature between both regenerator sides so for this maximal 
parameters we obtain a low difference of temperature. 
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Table 7 

Optimization criteria of the two responses. 



Difference of temperature 
between both regenerator 
sides (°C) 

Mechanical brake power 
for the Gamma Stirling 
engine (W) 

Objectif 

Minimise 

Maximise 

Minimal value 

20 

200 

Maximal value 

25 

290 

Cible 

20 

300 

Importance 

1 

2 


With consideration of the optimization criteria summarized in 
Table 7, the optimal functioning parameters which give the maximal 
brake power (290 W) and the minimal difference of temperature 
between both regenerator sides (20 °C) are the following: 

- Cooling water flow rate: 8.1 1/min. 

- Initial filling pressure: 8 bar. 

- Heating temperature: 350 °C. 

- Operation time: 7 min. 

A comparison with the results summarized in Table 6, show 
that the multi-objectif optimization has changed the initial optimal 
parameters ( ^able 6) in order to satisfy all imposed criteria ( Table 7). 

7. Concluding remarks 

In this paper, a new phenomenon causing the thermal energy 
dissipation in a Gamma type Stirling engine was studied. The 
asymmetry of temperature between both regenerator sides 
consumes a part of the exchanged thermal energy. The non- 
symmetrical design of Gamma Stirling engine is the principal 
causes of this phenomenon. An experimental design was per¬ 
formed to determine the influence of some operation parameters 
on the observed phenomenon. The heating temperature is the most 
significant parameter. An increase of 200 °C of the heating temper¬ 
ature increases the difference of temperature between both regen¬ 
erator sides of 20 °C. An analytical model linking the asymmetry of 
temperature to influent factors and interactions was proposed. 
With modeling, optimal operation parameters can be established. 
A minimal difference of temperature between both regenerator 
sides of 17 °C can be obtained by setting a minimal heating tem¬ 
perature and filling pressure and maximal cooling water flow rate. 

Second experimental design was proposed for a multi- objective 
optimization. The choosing Criteria are to maximize the engine 
power and to minimize the difference of temperature between 
both regenerator sides. A maximal power of 260 W and a differ¬ 
ence temperature of 20 °C are the optimal obtained responses. 

A geometric solution can be proposed to minimize the 
dissymmetry of temperature between both regenerator sides. The 
alimentation of engine cooler by the cooling water is actually made 
in only one side. A new cooler will be manufactured with alimen¬ 
tation in his both sides. The new configuration will ensure a homo¬ 
geneous distribution of working fluid temperature at input of the 
regenerator. 
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